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By means of oxide molecular beam epitaxy with shutter-growth mode, we have fabricated a series of electron- 
doped (Sri_;(;La;c) 2 lr 04 (001) (x = 0, 0.05, 0.1 and 0.15) single crystalline thin films and then investigated the 
doping dependence of electronic structure utilizing in-situ angle-resolved photoemission spectroscopy. We find 
that with increasing doping proportion, the Fermi levels of samples progressively shift upward. Prominently, an 
extra electron pocket crossing the Fermi level around the M point has been evidently observed in 15% nominal 
doping sample. Moreover, bulk-sensitive transport measurements confirm that doping effectively suppresses 
the insulating state with respect to the as-grown Sr 2 lr 04 , though doped samples still remain insulating at low 
temperatures due to the localization effect possibly stemming from disorders including oxygen deficiencies. Our 
work provides another feasible doping method to tune electronic structure of Sr 2 lr 04 . 

PACS numbers: 71.20.-b, 71.30.+h, 73.20.r, 77.55.Px 


5d transition-metal oxides (TMOs) have recently drawn a 
lot of attention. In these compounds, the complicated inter¬ 
play between spin-orbit coupling (SOC) and electron corre¬ 
lations has been suggested to host multiple novel quantum 
states, including topol ogic al Mott insulatorWeyl semimet- 
alJ^, axion insulatorJ^, and spin liquidJ^. These stud¬ 
ies were mostly initialized by the pioneering experiments on 
the prototype layered perovskite Sr2lr04, in which the strong 
SOC was revealed to lift the orbital degeneracy and then re¬ 
sult in a narrow half-filled Jeff = 1/2 band that even the rela¬ 
tively weak Coulomb repulsion of 5d electrons cou ld induce 
a Mott metal-insulator transition (MIT) thereiiPl^. Scatter¬ 
ing experiments discovered that, this Mott insulator, show¬ 
ing an effective pseudospin 1/2 antiferromagnetic (AFM) or¬ 
der at low temperature, can be well described by the Heisen¬ 
berg model with an exchange coupling of 60 to 100 meVll^. 
Such findings indicate that the low-energy behavior of this 
single-layer iridate rather resembles that of cuprates. Con¬ 
sequently, it is tempting to investigate the carrier doping of 
Sr2lr04, which might pave the way to discover a new family 
of unconventional superconductors. Indeed, a recent theoret¬ 
ical work has predicted the unconventional superconductivity 
in the electron-doped Sr2lrOi^. 

In this context, various experimental attempts have 
been performed to achieve the effective carrier doping of 
Sr2lrOzP^M^. So far, one of the most promising progresses is 
the try of surface electron doping of Sr2lr04 via in-situ potas¬ 
sium deposition, in which Fermi arcs and pseudogap behav¬ 
ior have been reported^, nevertheless, this finding is not con¬ 
sistent with the reports on the sibling La-doped Sr3lr207, in 
which only s mall Fermi pockets rather than Fermi arcs were 
reveale(p223. This controversy naturally raises one open ques¬ 
tion, namely, whether the existence of Fermi arcs is an univer¬ 
sal property of the electron-doped iridates irrespective of the 
specific doping method. As there has been no equivalent study 
reported so far for La-doped Sr2lr04, this important question 
still remains unanswered. 


In this letter, we synthesize a series of high-quality La- 
doped Sr2lr04 thin films on SrTiO3(001) single crystal sub¬ 
strates by means of oxide molecular beam epitaxy (0MBE), 
by which the metastable phases with large electron doping 
levels could be stabilized in the formation of thin films. We 
can thus investigate the electron-doping dependence of elec¬ 
tronic structure of bulk Sr2lr04 through in-situ angle-resolved 
photoemission spectroscopy (ARPES). We find that, with in¬ 
creasing the La-doped proportion, the Eermi levels of samples 
gradually shift upward to the conduction band, which is con¬ 
sistent with the electron doping. In addition, we recognize an 
extra electron-like band crossing the Eermi level upon 15% 
nominal La-doping, forming one pocket-like feature around 
the boundary of the reduced Brillouin zone (BZ), in stark con¬ 
trast to as-grown Sr2lr04 samples. Although electric transport 
measurements manifest the lanthanum substitution of Sr can 
significantly suppress the resistivity of this system, the doped 
Sr2lr04 still show clear insulating behavior. We argue that 
this may be caused by the random carrier hoping between lo¬ 
calized states in our system due to inevitable oxygen related 
defects. 

The La-doped Sr2lrO4(001) films were deposited onto 
SrTiO3(001) single crystalline substrates using a DC A R450 
OMBE system. The flat stepped surface of substrates was ob¬ 
tained by etching in hydrofluoric acid buffer solution (BHE) 
and annealing in a muffle furnace according the standard pre- 
cedurJ^. The in-plane lattice constant of SrTi03 substrates is 
3.905 A, which leads to a 0.44% tensile strain to the pseudo- 
cubic (La,Sr)2lr04 (3.888 A). The shuttered growth mode was 
applied to synthesize films in the distilled ozone atmosphere 
of 2 X 10 ^ Torr. During the growth, the temperature of 
substrates was kept at 800 °C verified by the thermocouple 
behind the sample stage. Moreover, the overall growth rate 
and surface structure of thin films were monitored by in-situ 
reflection high-energy electron diffraction (RHEED) during 
growth. Strontium (lanthanum) and iridium were evaporated 
from Knudsen effusion cells and an electron beam evaporator, 
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FIG. 1: (Color online) (a) (c) RHEED patterns of SrTiOsCOOl) sub¬ 
strate and as-grown (Sro. 95 Lao.o 5 ) 2 lr 04 along [110]^ azimuth direc¬ 
tion, respectively, (b) RHEED intensity curves integrated within 
the rectangle windows in panels (a) and (c). (d) Time depen¬ 

dence of the intensity of [00] diffraction streak in RHEED pattern 
of (Sri_;(;La;(;) 2 lr 04 (x = 0.05) thin film during deposition, (e) AEM 
topographic image of as-grown (Sro. 95 Lao.o 5 ) 2 lr 04 thin film, (f) The 
in-plane length dependence of height along the arrow direction in 
panel (e). 


respectively. All the doping levels as mentioned in this article 
are nominal values. 

Figures (a) and (c) show the RHEED patterns of 
SrTiO 3 ( 001 ) substrate before the growth and after the growth 
of 6 unit-cell (Sro. 95 Lao.o 5 ) 2 li '04 thin film, respectively, which 
were both taken with a glancing electron beam parallel 
to the [110]^ azimuthal direction. These RHEED images 
exhibit prominent Kikuchi lines, indicating the high crys¬ 
talline perfection of films. We compared the RHEED inten¬ 
sity curves integrated within the rectangle windows of the 
(Sro. 95 Lao.o 5 ) 2 lt '04 film and substrate [Eig.[ 2 b)], and the neg¬ 
ligible shift of the intensity peak indicates there is no evident 
lattice relaxation for the strained film. The amplitudes and 
periods of RHEED intensity oscillations of the [00] diffrac¬ 
tion streak tend to be constant once the Sr(La)/Ir flux ratio is 
adjusted to be close to stoichiometric 2. However, to further 
confirm that the shuttered SrO/(La 2 03 ) 0.5 and Ir 02 monolay¬ 
ers are complete during the growth, we measured the total 
thickness of Aims and their c-axis lattice constants through 
X-ray diffraction (XRD) to estimate the actual number of unit 
cells. By comparing this value to the number of shutting pe¬ 
riods of our calibration sample, we can then deduce the scal¬ 
ing factor, by which we can calibrate the Sr/La and Ir shut¬ 
ter opening times. In this way, the layer-by-layer growth of 
(Sri_;cLajc) 2 lr 04 Aims could be achieved. 

Eigure ^d) shows a typical RHEED intensity oscillation 
of the [ 00 ] diffraction streak (surrounded by a yellow pane 
in Eig. [^c)) of such a calibrated film. By means of the 
ex-situ atomic force microscopy (AEM), we further exam¬ 
ined the typical surface morphology of our thin film, as dis¬ 
played in Eig. [^e). The film surface shows a well defined 
terrace-step structure with a typical step height of 0.394 nm 


(a) Sr/La • Ir . Ti (b) (c) 




FIG. 2: (Color online) (a), (b) TEM cross-sections of Sr 2 lr 04 and 
(Sro. 85 Lao.i 5 ) 2 lr 04 thin films on SrTiO 3 ( 001 ) substrates, respectively, 
(c) Rocking curves of (006) diffraction peaks of Sr 2 lr 04 (orange 
symbol-line) and (Sro. 85 Lao.i 5 ) 2 lr 04 (blue symbol-line) thin films, 
respectively. Green lines are Gauss curves for fitting the rocking 
curves of thin films, (d) X-ray diffraction 6-26 scans of 8-nm-thick 
(Sri_;cLa;(:) 2 lr 04 (x = 0 and 0.15) thins films on SrTiO3(001) sub¬ 
strates. The characteristic peaks of substrates are labeled with an 
asterisk. 


[Fig.[ 2 f)], preserving the flat terraces and surface morphology 
of SrTiO 3 ( 001 ) substrates. This finding confirms the layer-by- 
layer growth mode of our Aims. Moreover, since all growths 
of our Aims were terminated after an integral number of oscil¬ 
lation finishing, we can expect a single Ir 02 surface termina¬ 
tion for all the Aims. 

To characterize the phase purity and crystallinity of 
(Sri_;cLa;c) 2 lr 04 thin Alms, the high-resolution cross-sectional 
transmission electron microscopys (TEM) and the ex-situ 
XRD measurements were performed. Eigs. 2(a) and (b) show 
the TEM images of the typical Sr 2 lr 04 and (Sro. 85 Lao.i 5 ) 2 E 04 
thin Alms deposited on SrTiO 3 ( 001 ) substrates. Both data 
demonstrate the continuous and evenly spaced SrO/(La 2 03 ) 0.5 
and Ir 02 layers and the atomically sharp interfaces between 
the epitaxial thin Alms and the substrates, with no discernible 
interdiffusion or vacancy observable. The XRD results of typ¬ 
ical (Sri_;cLa;c) 2 lr 04 thin Alms are shown in Eigs. |^c) and 
(d). The 6-26 scans of the Alms [Eig.j^d)] are consistent with 
the growth of phase-pure (OOl)-oriented (Sri_;cLa;c) 2 li' 04 - The 
clear Kiessig interface fringes indicate smooth Aim surfaces 
and sharp interfaces between the thin Alms and substrates. We 
also measured the rocking curves around the (006) peaks of 
Sr 2 lr 04 and (Sro. 85 Lao.i 5 ) 2 E 04 , as shown in Eig. [^c). The 
full widths at half maximum (EWHM) of these curves are 
0.062° and 0.069°, respectively, indicating the good epitax¬ 
ial Aim quality. In addition, through Bragg’s function on the 
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FIG. 3: (Color online) (a) Photoemission intensity plots of Sr 2 lr 04 thin film along X-T-X direction, taken with 21.2 eV photon energy at 40 
K. (b) The valence-band photoemission spectra at F point for (Sri_;(-La;(;) 2 lr 04 thins films, (c) Photoemission intensity plots along X-M of 
Sr 2 lr 04 thin film, (e), (g) Photoemission intensity plots along X-M and M-T of (Sro.85Lao.i5)2li'04, respectively, (d), (f), (h) The momentum 
distribution curves for the data in panels (c) , (e) and (g). (i) The two-dimensional projection of the Brillouin zone and the high-symmetry 
directions. The blue solid lines represent the reduced Brillouin zone, (j) Fermi surface map integrated over [Ef - 15 meV, Ep -i- 15 meV] of 
(Sro.85Lao.i5)2lr04 thin film. Yellow dashed lines guide for the elliptical Fermi surface around M. BLue dashed lines illustrate the refiected 
shadow band of Fermi arcs. All the APRES data were taken with 21.2 eV photon energy at 40 K. 


(006) peak positions, we can deduce the out-of-plane c lattice 
constants of Sr 2 lr 04 and (Sro.85Lao.i5)2lr04 to be 12.782A and 
12.821 A, respectively. For comparison, the out-of-plane c lat¬ 
tice constant of bulk Sr 2 lr 04 was determined to be 12.899Ap^. 
Obviously, the in-plane tensile strain for the epitaxial films re¬ 
sults in the compression in the out-of-plane c direction, while 
the doping of lanthanum with smaller ionic radius partially 
counterbalances this effect. From another perspective, this 
finding confirms the real bulk doping of lanthanum in our 
(Sri_;cLajc) 2 lr 04 films grown by OMBE. 

After growth, we performed in-situ ARPES measurements 
on these samples in order to investigate the electron doping 
evolution of (Sri_;cLa;c) 2 lr 04 - Thin films were transferred 
through an ultrahigh vacuum buffer chamber (1.0 x 10“^^ torr) 
to the combined ARPES chamber for measurements immedi¬ 
ately after the growth. This ARPES system is equipped with 
a VG-Scienta R8000 electron analyzer and a SPECS LIVES 
helium discharging lamp. The data were collected at 40 K un¬ 
der ultrahigh vacuum of 8 x 10“^^ torr. The angular resolution 
was 0.3 and the overall energy resolution was set to 15 meV 
(Hel, 21.2 eV photon energy). During the measurements, the 
films were stable and did not show any sign of degradation. 

Eig. I^a) shows the valence band photoemission intensity 
plot for Sr 2 lr 04 epitaxial films along X-T-X high-symmetry 
direction highlighted by the red line in Eig. 3(i). By compar¬ 
ing with DET calculations, we can identify that the features 
between - 7.0 and - 2.0 eV are mainly contributed by the O 
2p states, while the Ir t 2 g orbitals mainly distribute the feature 
located at around - 0.5 eV, as illustrated b y the integrated spec¬ 
tra in momenta around E point [Eig.[^b)P^E3f^. While, upon 
further La doping (from x = 0 to 0.15 nominal doping level), 
we discovered a continuous energy shift (as large as 500 meV) 
for these spectral features as shown in Eig. |^b). This find¬ 


ing confirms the gradual electron doping into the insulating 
Sr 2 lr 04 by substituting Sr with La, which would increase the 
chemical potential of this system and keep pushing valence 
band (VB) features to higher binding energy. Here, all spectra 
have been renormalized by the spectral intensity at the binding 
energy of 8.0 eV for comparison. 

Surprisingly, as for the sample with nominal 15% La- 
doping, we observe an extra electron-like band crossing 
the Eermi level [Eigs. Se) and (g)], while the non-doped 
Sr 2 lr 04 just presents as a typical insulator with an energy 
gap[Eig. I^c)]. Moreover, the comparison between the cor¬ 
responding momentum distribution curves in Eigs. [^d), (f) 
and (h) further confirms this finding. As shown in Eig. [^), 
the band forms one elliptical electron pocket around the M 
point of the reduced Brillouin zone which is guided by yel¬ 
low dashed line[Eig. [^i)]. It seems that our findings are in 
qualitative agreement with the Eermi pockets discovered in 
the metallic electron-doped Sr 3 lr 207 , in which electron carri¬ 
ers are doped into the conduction band and produces a small 
Eermi surface containing only the added x electronJ^. How¬ 
ever, we could still not simply rule out the possibility of the 
Eermi arcs. Taking into account of the significant octahedral 
rotation and the possible antiferromagnetic (AE) ordering in 
(Sri_;cLa;c) 2 li' 04 , it is highly possible for ARPES to observe 
the reflected shadow bands of the Eermi arcs, which would as 
well produce the pocket-like feature as illustrated by the blue 
dashed lines in Eig.[^). To further pin down the origin of this 
pocket-like feature, a more detailed doping dependence of the 
low-lying electronic structure is highly imperative. 

To further determine whether our findings are caused by 
the surface or the bulk electronic structure, we compared the 
in-plane resistivity of Aims with different doping levels. As 
shown in Eig.Qa), for the as-grown Sr 2 lr 04 , the normalized 
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FIG. 4: (Color online)(a) The temperature dependence of resis¬ 
tance for (Sri_;(;.La;c)2lrO4(001) (x = 0 and 0.15) epitaxial films, 
normalized by the resistance at 300 K temperature, (b) The loga¬ 
rithm of the in-plane normalized resistivity as a function of 
for (Sro.85,Lao.i5)2lr04 epitaxial film. Inset: Sketch of a typical two- 
dimensional variable range hoping model. 

R-T curve manifests a typical insulator behavior and the re¬ 
sistivity has been out of measurement range when the temper¬ 
ature is below 70 K; nevertheless, though the film with 15% 
nominal doping level still behaves as an insulator, its in-plane 
resistivity has dropped by around three orders at 70 K, and it 
is still in the range of our measuring system at as low as 5 
K. This prominent change of the bulk properties of films is in 
good agreement with our ARPES result, and it undoubtedly 
supports that our findings are mainly caused by the electron 
doping into the bulk Sr 2 lr 04 . We note that (Sro.85Lao.i5)2lr04 
is still insulating though there is finite spectra weight in the 
vicinity of Fermi energy (Ep), which implies that this as¬ 
tern might be in the strongly localized regimJ^. Fig. Qb) 
shows the logarithm of the in-plane resistivity as a function 
of T“^/^ for a temperature range from 5 to 300 K. In the low 
temperature range (5 to 40 K), we could observe an excellent 
agreement with a linear fit, consistent with the typical two¬ 


dimensional variable range hopping (VRH) modei^. In this 
case, carriers would hop between localized states as sketched 
by the inset of Fig. [^b). At higher temperatures, though 
there exists a slight deviation from the linear fit, we can con¬ 
clude that the VRH may still be the main conduction mecha¬ 
nism. Recently, one scanning tunneling microscopy and spec¬ 
troscopy experiment discovered some oxygen related defects 
on the surface of Sr 2 lr 04 , which are likelymative to the sample 
and regardless of the good sample qualitySl. Such random dis¬ 
tribution of the defects are in good agreement with the VRH 
behavior in our in-plane resistivity data. 

In summary, we have grown a series of La-doped 
(Sri_;cLa;c)2lrO4(001) (x = 0, 0.05, 0.1, and 0.15) epitaxial 
films on SrTiO3(001) substrates and studied the correspond¬ 
ing electronic structure and transport properties. By increas¬ 
ing the La-doping level, Fermi level for the samples are drifted 
up getting close to higher binding energy. Moreover, a fast 
dispersing state is observed around M point for the sam¬ 
ples with 15% La-doping. In addition, the R-T curves of 
(Sro. 85 Lao.i 5 ) 2 F 04 thin film is suppressed under the one of 
undoped Sr 2 lr 04 thin film. Although the La-doped propor¬ 
tion is as high as 15%, the films still remain insulators, which 
might be resulted from the strong localized effect induced by 
disorders including the oxygen defects in the samples. 

We gratefully acknowledge the helpful discussion with 
Prof. Donglai Feng, Dr. Rui Peng and Dr. Wei Li. This work 
was supported by National Basic Research Program of China 
(973 Program) under the grant Nos. 2011CBA00106 and 
2012CB927400, the National Science Foundation of China 
under Grant Nos. 11274332 and 11227902, and Helmholtz 
Association through the Virtual Institute for Topological In¬ 
sulators (VITI). M. Y. Li and D. W. Shen are also supported 
by the ’’Strategic Priority Research Program (B)” of the Chi¬ 
nese Academy of Sciences (Grant No. XDB04040300). 


* Electronic address: dwshen@mail.sim.ac.cn 
^ D. Pesin and L. Balents, Nat. Phys. 6, 376 (2010). 

^ A. A. Burkov and L. Balents, Phys. Rev. Lett. 107,127205 (2011). 
3 K. Yu. Yang, Y. M. Lu, and Y. Ran, Phys. Rev. B 84, 075129 
( 2011 ). 

W. Witczak-Krempa, G. Chen, Y. B. Kim, and L. Balents, Ann. 
Rev. Cond. Matt. Phys. 5, 57-82 (2014). 

^ X. G. Wan, A. Vishwanath, and S. Y. Savrasov, Phys. Rev. Lett. 

108, 146601 (2012). 

^ L. Balents, Nature 464, 199 (2010). 

^ X. G. Wen, Phys. Rev. B 65, 165113 (2002). 

^ M. J. Lawler, A. Paramekanti, Y. B. Kim, and Leon Balents, Phys. 
Rev. Lett. 101, 197202 (2008). 

^ B. J. Kim, H. Ohsumi, T. Komesu, S. Sakai, T. Morita, H. Takagi 
and T. Arima, Science 323, 1329 (2009). 

B. J. Kim, Hosub Jin, S. J. Moon, J.-Y. Kim, B.-G. Park, C. S. 
Leem, Jaejun Yu, T. W. Noh, C. Kim, S.-J. Oh, J.-H. Park, V. 
Durairaj, G. Cao, and E. Rotenberg, Phys. Rev. Lett. 101, 076402 
(2008). 

11 Jungho Kim, D. Casa, M. H. Upton, T. Gog, Young-June Kim, J. 
F. Mitchell, M. van Veenendaal, M. Daghofer, J. van den Brink, 


G. Khaliullin, andB. J. Kim, Phys. Rev. Lett. 108, 177003 (2012). 
1^ S. Fujiyama, H. Ohsumi, T. Komesu, J. Matsuno, B. J. Kim, M. 
Takata, T. Arima, and H. Takagi, Phys. Rev. Lett. 108, 247212 
( 2012 ). 

1^ H. Watanabe, T. Shirakawa, and S. Yunoki, Phys. Rev. Lett. 110, 
027002 (2013). 

i"i M. Ge, T. F. Qi, O. B. Korneta, D. E. De Long, P. Schlottmann, 
W. P. Crummett, and G. Cao, Phys. Rev. B 84, 100402 (2011). 

1^ T. F. Qi, O. B. Korneta, L. Li, K. Butrouna, V. S. Cao, Xiangang 
Wan, P. Schlottmann, R. K. Kaul, and G. Cao, Phys. Rev. B 86, 
125105 (2012). 

1^ J. Nichols, J. Terzic, E. G. Bittle, O. B. Korneta, L. E. De Long, J. 
W. Brill, G. Cao, and S. S. A. Seo, Appl. Phys. Lett. 102, 141908 
(2013). 

1^ C. R. Serrao, J. Liu, J. T. Heron, G. Singh-Bhalla, A. Yadav, S. J. 
Suresha, R. J. Pauli, D. Yi, J.-H. Chu, M. Trassin, A. Vishwanath, 
E. Arenholz, C. Frontera, J. Zelezny, T. Jungwirth, X. Marti, and 
R. Ramesh, Phys. Rev. B 87, 085121 (2013). 

L. D. Miao, H. Xu, and Z. Q. Mao, Phys. Rev. B 89, 035109 
(2014). 

J. S. Lee, Y. Krockenberger, K. S. Takahashi, M. Kawasaki, and 











5 


Y. Tokura, Phys. Rev. B 85, 035101 (2012). 

O. B. Korneta, T. F. Qi, S. Chikara, S. Parkin, L. E. De Long, P. 
Schlottmann, and G. Cao, Phys. Rev. B 82, 115117 (2010). 

Y. K. Kim, O. Krupin, J. D. Denlinger, A. Bostwick, E. Rotenberg, 
Q. Zhao, J. E. Mitchell, J. W. Allen, and B. J. Kim, Science 345, 
6193 (2014). 

A. de la Torre, E. C. Hunter, A. Subedi, S. McKeown Walker, A. 
Tamai, T. K. Kim, M. Hoesch, R. S. Perry, A. Georges, and E. 
Baumberger, Phys. Rev. Lett. 113 , 256402 (2014). 

J.H. He, H. Hafiz, T. R. Mion, T. Hogan, C. Dhital, X. Chen, Q. 
Lin, M. Hashimoto, D. H. Lu, Y. Zhang, R. S. Markiewicz, A. 
Bansil, S. D. Wilson, R.-H. He, Scientific Reports 5, 8533 (2015). 
2"^ M. Kawasaki, K. Takahashi, T. Maeda, R. Tsuchiya, M. Shino- 
hara, O. Ishiyama, T. Yonezawa, M. Yoshimoto, and H. Koinuma, 
Science 266, 1540(1994). 


23 M. K. Crawford, M. A. Subramanian, R. L. Harlow, J. A. 
Eernandez-Baca, Z. R. Wang, and D. C. Johnston, Phys. Rev. B 
49 , 9198 (1994). 

2^ R. Arita, J. Kunes, A. V. Kozhevnikov, A. G. Eguiluz, and M. 

Imada, Phys. Rev. Lett. 108 , 086403 (2012). 

22 C. Martins, M. Aichihorn, L. Vaugier, and S. Biermann, Phys. 
Rev. Lett. 107 , 266404 (2011). 

28 H. Watanabe, T. Shirakawa, and S. Yunoki, 105, 216410 (2010). 
2^ R. Scherwitzl, S. Gariglio, M. Gabay, P. Zubko, M. Gibert, J.-M. 

Triscone, Phys. Rev. Lett. 106 , 246403 (2011). 

3^ W. Brenig, G. H. Dohler, and H. Heyszenau, Philos. Mag. 27, 
1093 (1973). 

31 J.X. Dai , E. Calleja , G. Cao, K. McElroy, Phy. Rev. B 90 , 
041102(R) (2014). 



